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Abstract 

Based on CORSIKA and REAS2 simulations, we investigate the dependence of 
geosynchrotron radio emission from extensive air showers on the energy of the pri- 
mary cosmic ray and the depth of the shower maximum. It is found that at a 
characteristic lateral distance, the amplitude of the bandpass-filtered radio signal 
is directly proportional to the energy deposited in the atmosphere by the electro- 
magnetic cascade, with an RMS uncertainty due to shower-to-shower fluctuations 
of less than 3%. In addition, the ratio of this radio amplitude and that at a larger 
lateral distance is directly related to the atmospheric depth of the shower maxi- 
mum, with an RMS uncertainty of ~ 15-20 gem -2 . By measuring these quantities, 
geosynchrotron radio emission from cosmic ray air showers can be used to infer 
the energy of the primary particle and the depth of the air shower maximum on a 
shower-to-shower basis. 
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1 Introduction 



During the last few years, the technique of radio detection of cosmic ray air 
showers has experienced an impressive renaissance [H; 0; 0; 0] • The activities 
are driven by the prospect of establishing a new observing technique with 
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nearly 100% duty cycle and very good angular resolution, which would com- 
plement the existing particle detectors and air fluorescence telescopes. Due 
to the coherent nature of the radio emission, i.e. an approximately quadratic 
scaling of the emitted radio power with the energy of the primary cosmic ray, 
the technique is particularly well-suited for the detection of ultra-high energy 
cosmic rays. 

One important question is how the observables measured with the radio tech- 
nique can be related to the energy and mass of the primary cosmic ray particle. 
In this article, we use the geosynchrotron model for radio emission from cosmic 
ray air showers 0; @; Q; [i[ to demonstrate how shower parameters could be 
derived from radio-only measurements that allow one to estimate the primary 
energy and composition. It is shown that the characteristics of the lateral 
distribution function of the radio emission can be exploited to estimate on a 
shower-to-shower basis both the energy deposited in the atmosphere and the 
depth of the shower maximum. The energy determination can be performed 
with a measurement at an optimum lateral distance which is independent of 
the energy of the primary particle in the considered energy range from 10 18 
to 10 20 eV — it only depends on the air shower geometry and radio observing 
frequency. Combining this measurement with an additional measurement at a 
different lateral distance provides a handle on the depth of the shower maxi- 
mum and consequently yields information related to the mass of the primary 
particle. 



All results presented in this simulation study were derived with the geosyn- 
chrotron model as implemented in the REAS2 code (§]. Additional mechanisms 
which can contribute to the radio signal, such as Cherenkov-like emission from 
a charge excess and contributions from the net charge variation during the air 
shower evolution j§; flo| . are the subject of current investigations and are not 
included in the analysis presented here. A direct comparison between the im- 
plementations of different emission models by various authors is difficult due 
to the very different calculational approaches. In contrast to other available 
models, however, the geosynchrotron model implemented in the REAS2 code 
allows us an absolute, parameter-free calculation of the emission from indi- 
vidual air showers, accounting for realistic shower-to-shower fluctuations and 
detailed spatial, angular and energy distributions of the shower particles. If 
radio emission from air showers is dominated by geomagnetic effects — an 
assumption supported by experimental data showing strong correlations be- 
tween the strength of the radio signal and the orientation of the air shower 



axis to the Earth's magnetic field [1|; [III; ll2( — the results of this analysis can 
be considered generally applicable. 



In the following we first describe the methodology applied throughout the 
present analysis. We then investigate the case of air showers with 60° zenith 
angle, where the effects are very prominent and could be well exploited ex- 
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perimentally. Afterwards we demonstrate that the same qualitative behaviour 
is fulfilled for air showers with 45° zenith angle, where it would, however, be 
harder to use experimentally, followed by our conclusions. 



2 Methodology 



We investigate four air shower geometries: showers with 60° zenith angle com- 
ing from the south (180+90 showers), from the east (180 showers) and from 
the north (180 showers) as well as air showers with 45° zenith angle coming 
from the south (180 showers). Each set of 180 showers covers the energy range 
from 10 18 eV to 10 20 eV according to a power-law with spectral index of —1 and 
consists of 60 proton-induced, 60 iron-induced and 60 photon-induced show- 
ers. The additional set of 90 showers has been simulated for a fixed energy of 
10 19 eV, induced in equal parts by proton, iron and photon primaries. 



For each of these 810 showers, a complete simulation with CORSIKA 6.710 
using the hadronic interaction models QGS JETII-03 1% |l5j and UrQMD 
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1.3.1 |16|; Uj\ was performed to produce detailed histogram information to be 
used for the radio simulations as described in . The magnetic field was set to 
a strength of 0.23 Gauss with an inclination angle of -37.3°, corresponding to 
the values valid for the southern site of the Pierre Auger Observatory 18] in 
Argentina. Particle thinning within CORSIKA was set to 10 -6 with optimised 



weight limitation according to [19j]. Electromagnetic particles were tracked 



down to kinetic energies of 400 keV using EGS4 [20j. For photon-primaries 
with energies > 10 19 eV, pre-showering in the Earth's magnetic field was taken 
into account. Also, the LPM effect | 
is accounted for in the simulations. 



22|; [23j as implemented in CORSIKA 



In a second step, the radio emission for each of the 810 CORSIKA-simulated 
air showers was calculated with REAS 2.58 @;Q for a regular grid of observer 
locations at a height of 1400 m above sea level, again adequate for the southern 
Pierre Auger Observatory. For the 45° and 60° zenith angle showers coming 
from the south, in addition, the radio signal at sea level was calculated. 



The resulting unlimited bandwidth time-series pulses were then digitally fil- 
tered with three idealised rectangle filters: 16-32 MHz, 32-64 MHz and 64- 
128 MHz. (The filtered pulses correspond to signals that experiments with the 
given frequency bandwidths would measure; the bandwidth of 32-64 MHz is 
closest to the one employed in current experimental approaches.) The abso- 
lute values (vectorial sum of all three polarisation components) of the peak 
amplitudes of these filtered pulses, derived for each observer location and air 
shower, form the basis for the following discussion. 
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3 The 60° zenith angle case 



We first discuss the important characteristics of the radio signal for the case 
of 60° zenith angle. In a first step we shall consider showers of 10 19 eV energy 
and later show that the results also apply to the energy range from 10 18 to 
10 20 eV. 



3. 1 Characteristics of lateral profiles 



Earlier analyses [3] have shown that the slope of the radio lateral distribution 
function is correlated with the depth of the air shower maximum, X ma xj_and 



consequently contains information on the mass of the primary particle 25 
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Fig. 1. Lateral distribution of the 32-64 MHz filtered peak radio amplitude for 
10 19 eV showers coming from the south and observers north of the shower core. 
Estimates of peak radio amplitudes that would yield a signal-to-noise ratio (SNR) 
of 10 for ideal (galactic plus atmospheric), rural and urban noise as given in [26J 
are marked. Please note that the SNR scales with the filtered peak radio amplitude 
squared. 

In Fig. 1, we show the lateral distribution of the 32-64 MHz filtered peak 
radio amplitude derived for the 90 air showers with fixed energy of 10 19 eV. 
The lateral distance is given in ground coordinates^ in the direction defined 
by the continuation of the air shower axis (i.e., in this case to the north, as the 
showers are coming from the south). Measurements at distances up to one km 
should be feasible for ultra-high energy cosmic rays, as the comparison with 
continuous noise estimates based on ITU/CCIR measurements taken from 26 
illustrates. 



1 We use ground-distances rather than shower-distances throughout this article, as 
shower-distances would not remove the intrinsic asymmetries of the radio signal and 
ground distances are experimentally relevant in the end. 
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There is an important feature in these lateral distributions: while the slopes 
of the individual lateral profiles differ significantly, all profiles intersect in a 
small region around ~ 250 m. This essentially arises from the strong forward 
collimation of the radio emission in combination with the integration over the 
course of the air shower evolution. 
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Fig. 2. Peak filtered radio field 
strengths of 10 19 eV showers coming 
from the south as measured in the flat 
region 275 m north of the shower core. 
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Fig. 3. Same as Fig. 2, but in a steep 
region 725 m north of the shower core. 



A more quantitative view of the 32-64 MHz filtered peak amplitudes in this 
region and their relation to the individual shower X mscK values is given in Fig. 
2. The distribution of filtered field strengths is approximately independent of 
X miB — and thus almost independent of shower to shower fluctuations. We call 
the distance region of the intersection, characterised by the minimum in the 
RMS spread of the amplitudes, flat region from now on. (The location of this 
region depends on several parameters, in particular the shower zenith angle 
and window of observing frequency, but also the chosen normalisation applied 
to the radio amplitudes, as will be discussed in the following.) Although the 
amplitude distribution is approximately flat, there are visible "steps" between 
the values for the different types of primary particles. The reason for these 
"steps" is the fact that the radio emission is generated only by the electro- 
magnetic cascade of an air shower, and depending on the type of the primary 
particle, a different fraction of the energy of the primary particle is transferred 



to the electromagnetic cascade [27|; [28 



At large lateral distances, e.g. around 700 m, the lateral profiles for the indi- 
vidual air showers plotted in Fig. 1 lead to very different filtered electric field 
strengths. Such a region will be called steep region from now on. A quantita- 
tive look at the filtered electric field strengths as a function of A max is given in 
Fig. 3. In a steep region, the filtered electric field strength is directly correlated 
with the X max of the air shower. 
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3.2 Energy determination 



With experimental measurements, the arrival direction and core position of an 
air shower can be deduced from radio-only measurements 29|. The energy of 
the primary particle, on the other hand, is a priori unknown. The correlations 
discussed in the previous section can, however, be generalized for a range of 
energies. 



In Fig. 4, the lateral radio emission profiles of the 180 simulated air showers 
spanning the energy range from 10 18 to 10 20 eV are shown after filtering to 
the observing bandwidth of 32-64 MHz. To make shower signals comparable 
we have scaled them with the energy that the electromagnetic cascade of 
the air shower deposits in the atmosphere. While the absolute electric field 
amplitudes of the 180 showers span three orders of magnitude, the normalised 
lateral distribution profiles exhibit the same qualitative behaviour as those 
of the monoenergetic showers depicted in Fig. 1. In particular, there exists 
again a well-defined region for all showers (and thus all energies), in which the 
normalised lateral profiles intersect. As will be discussed later, the location of 
the intersection region and the spread of the normalised field strengths at this 
distinctive lateral distance depends on the quantity with which the profiles 
are scaled. 
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Fig. 4. Lateral distribution of the 32-64 MHz filtered peak radio amplitude nor- 
malised with the energy deposited in the atmosphere by the electromagnetic cas- 
cade. Shown are 180 air showers spanning the energy range from 10 18 to 10 20 eV. 
The normalised field strength profiles show the same behaviour as those of the mo- 
noenergetic showers in Fig. 1. In particular, all normalised lateral profiles intersect 
in a narrow region around ~ 175 m lateral distance. 



Already in 301]. we have demonstrated that the radio signal correlates with 
the sum of the number of electrons and positrons in the air shower maximum, 
Amax- This correlation is illustrated in Fig. 5, where the filtered electric field 
strength per iV max in the appropriate flat region of all 180 showers with 60° 
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zenith angle coming from the south, spanning the energy range from 10 18 eV 
to 10 20 eV, is plotted against X max . Regardless of the energy and type of 
the primary particle and the value of X max , the filtered peak electric field 
strength per iV max yields approximately the same value. In particular, the 
"steps" between the different types of primary particles seen in Fig. 2 have 
been reduced considerably by the normalisation with the maximum number 
of particles in the electromagnetic cascade. 

Another quantity with which the radio signal could be normalised is the 
"calorimetric energy" in the air shower, E r ^ . (This is the quantity with which 
the fluorescence light correlates best, cf. 28|;|3l|.) The calorimetric energy of 
each simulated air shower is calculated from the total energy deposited in an 
infinitely thick atmosphere. As demonstrated in Fig. 6, for 60° zenith angle 
showers, the calorimetric energy E ca \ has an even better correlation with the 
radio signal in the appropriate flat region than A" max . The reason is that E ca \ 
fluctuates less from shower to shower than the value of N^^. 
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Fig. 5. At 225 m north of the shower 
core, the radio pulse height per iV max 
over energies of 10 18 to 10 20 eV and all 
three particle types is approximately 
constant, i.e., this is a flat region. 
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Fig. 6. Normalising the radio signal 
with the calorimetric energy E ca \ of an 
air shower, yields even better results. 
The flat region then lies at 175 m north 
of the shower core. 
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From a theoretical point of view, however, the quantity with which the radio 
signal should correlate best is the "energy deposited in the atmosphere by the 
electromagnetic cascade" , E em . In contrast to E ca x, E em neither includes energy 
deposited by muons or hadrons (which do not contribute to the radio signal) 
nor takes into account the energy deposited by particles in the ground. The 
high quality of the correlation with E em is illustrated in Fig. 7, which shows 
the normalised radio amplitudes at 175 m north of the shower core, where the 
corresponding flat region lies (see also Fig. 4). 

In Fig. 8 we compare the width of the correlation with the three discussed 
quantities, A" max , E CB x and E em . As mentioned earlier, the location of the flat 
region, i.e., the minimum of the RMS distribution, depends on the chosen 
normalisation. In addition, it becomes clear that E em yields the best correla- 
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Fig. 7. For radio pulse heights nor- 
malised with E em , the flat region lies 
again at a distance of ~ 175 m north of 
the shower core. 
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Fig. 8. The lowest deviation from the 
mean normalised amplitude is reached 
for the normalisation with E em with 
an RMS of < 3% in the flat region at 
~ 175 m north. 



tion, with an RMS deviation from the mean normalised amplitude of <3% at 
a ground distance of ~ 175 m. We therefore apply the normalisation of the 
radio signal with E em in the following. 
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Fig. 9. A measurement in the flat re- 
gion directly yields the shower E em on 
a shower-to-shower basis via a linear re- 
lation. The fit paramater is a = 197900. 
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Fig. 10. The RMS spread in peak field 
strengths per E em in the flat region over 
all energies and particle types amounts 
to <3%. 



Once the value of the normalised, energy-independent radio amplitude per E em 
in the flat region is known for a given geometry and observing frequency - 
either from simulations or from measurements — a measurement of the peak 
filtered radio amplitude y4 flat in the flat region can be used to immediately 



read off the value of E em in a linear relation, 



E em = a y4 fla t yuV 1 m GeV 



where a is a shower geometry dependent parameter, see Fig. 9. In Fig. 10, 
the distribution of the 32-64 MHz filtered peak radio amplitudes per E em 
around their mean value is illustrated. As discussed earlier, the RMS of this 
distribution amounts to <3%. In principle, i.e. not considering experimental 
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uncertainties, the E em value can therefore be inferred with a theoretical uncer- 
tainty of <3%. It should be stressed that this theoretical accuracy would be 
achievable on a shower-to-shower basis, unhindered by shower-to-shower fluc- 
tuations and without a priori knowledge of the mass of the primary particle. 
The derived value of E em can in turn be related to the ener gy of the primary 



particle with an additional uncertainty of the order ~5-7% 32 



3.3 Composition determination 



The ratio of the measured field strengths in the flat region and steep region 
can be used to infer the X max value on a shower-to-shower basis, which can in 
turn be related to the mass of the primary particle. This is demonstrated in 
Fig. 11. The correlation of X max with the ratio A^ t /A steep can be described 
well with the fit function 
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The deviations of the X max values of individual showers from this fit are de- 
picted in Fig. 12. The RMS spread of these deviations for the 180 showers 
shown here amounts to 15.9 gem -2 . 



Note that this spread only accounts for intrinsic shower-to-shower fluctuations 
and does not include any measurement errors. For comparison, fluorescence 
telescopes, which measure directly, reach an experimental resolution of 

the order of 35 to 20 g cm -2 [3c 
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Fig. 11. The ratio of peak field 
strengths in the flat region and steep 
region yields direct information on the 
X max value of an individual air shower. 
The fit parameters are a = 856.1, 
b = 0.3149 and c = 0.4340. 
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Fig. 12. The RMS spread of the indi- 
vidual shower X max values around the 
values calculated with the fit function 
amounts to 15.9 gem -2 . 
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3.4 Influence of observing frequency 



The characteristic lateral distance at which the flat region is located depends 
strongly on the zenith angle of the air shower. The frequency window of the 
radio observations, however, also influences its location. For the case of 60° 
zenith angle air showers and normalisation with E cm , this is illustrated in 
Fig. 13. A clear minimum of the RMS with values <3% can be deduced for 
each frequency band separately. The corresponding distance of the flat region 
amounts to the aforementioned ~ 175 m in case of the 32-64 MHz band. For 
higher frequencies, it lies at smaller distances, while for lower frequencies, it 
shifts to larger distances. From an experimental point of view, larger distances 
are favourable because they require less densely equipped radio antenna arrays. 
Multi-frequency observations including lower observing frequencies would thus 
allow an even better exploitation of measurements in the flat and steep regions 
for the determination of E em and X max . 
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Fig. 13. RMS spread of the peak field 
strengths per E em as a function of 
ground distance from the shower core 
in the azimuthal direction defined by 
the shower axis for observers at 1400 m 
above sea level. 
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Fig. 14. Same as Fig. 13 but for 
observers at sea level. The flat regions 
shift to slightly larger distances. 



3.5 Influence of observation altitude 



The position of the flat region is also influenced by the observer's height above 
sea level. This is due to the fact that the geometrical distance from the shower 
maximum to the radio antennas changes with the observer position, which 
influences the lateral distributions of the radio emission [3]. 

In Fig. 14, the RMS deviation from the mean filtered radio amplitude per E eva 
is shown for observers at sea level, in contrast to the 1400 m above sea level 
used as standard value in this analysis. While the difference to an observer at 
1400 m above sea level (Fig. 13) is not dramatic for 60° zenith angle showers, 
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there is a clear trend for the flat regions to move to larger lateral distances, 
as can be expected from the flattening of the lateral profiles. 

3.6 Influence of observer azimuth angle 

Radio emission from inclined air showers exhibits a significant asymmetry for 
observer positions along the azimuthal direction given by the air shower axis 
and positions perpendicular to the air shower axis. For inclined air showers, 
part of this asymmetry is caused by geometrical effects and can be taken out 
by changing from ground coordinates to shower coordinates [3]. The geosyn- 
chrotron model, however, alsopredicts an intrinsic asymmetry that is present 
even for vertical air showers [8j. In Fig. 15, we illustrate the RMS spread of 
the 32-64 MHz filtered peak amplitudes per E em for observers along different 
azimuthal directions in case of air showers coming from the south. The be- 
haviour is almost identical for observers along the two azimuthal directions 
defined by the shower axis, in this case to the north and to the south, except 
for a small asymmetry caused by the inclination of the geomagnetic field. 

For an observer in the north-west, i.e. at 45° azimuth to the shower axis, the 
flat region shifts to a very small lateral distance, which is no longer resolved 
in the simulations. This is the behaviour expected from the model-predicted 
asymmetries of the radio "footprint" . In addition, there is a local minimum in 
the RMS distribution around ~ 725 m. On a closer look, such additional local 
minima are also visible in the north direction, especially at 64-128 MHz (cf. 
Fig. 13). These local minima are dependent on the window of observing fre- 
quency and therefore related to interference effects arising from the bandpass 
filtering. 

For the east or west directions, i.e. perpendicular to the shower axis, the flat 
region also lies at distances too small to be resolved by the simulations. At 
the same time, additional prominent local minima occur at large distances. 

From an experimental point of view, large lateral distances for the flat re- 
gion are preferred, and measurements would thus be easiest in the azimuthal 
directions given by the air shower axis. 

3. 7 Influence of shower azimuth angle 

So far, we have only treated the case of 60° zenith angle showers coming 
from the south. The qualitative behaviour of the radio signal, however, is 
independent of the shower azimuth angle. This is illustrated in Fig. 16, where 
the peak filtered radio amplitude per E ein in the flat region is plotted for the 



11 



50 r- 
45 - 
40 - 
35 - 
30 - 
25 - 
20 - 
15 - 
10 

5 







north — h 

south ----x 

east * 
west 

northwest ■■■ ■■ 




100 200 300 400 500 600 700 800 900 1000 
distance from shower core [m] 

Fig. 15. RMS spread of the 32-64 MHz peak field strengths per E em for different 
observer directions relative to the shower core. The local minima are dependent on 
the window of observing frequency and thus related to the bandpass filtering. 

three sets of air showers coming from the south, east and north. (The flat region 
is then, correspondingly, found at ~ 175 m to the north, west and south of the 
shower core.) While the absolute value of the filtered peak radio amplitude per 
E em varies with the angle between the shower axis and the Earth's magnetic 
field (called geomagnetic angle, corresponding to 67.3°, 107.6° and 172.7° for 
the three cases), the values for a given geometry are again approximately equal 
over all showers. Once this constant value is known for a given geometry - 
from simulations or measurements — the E em determination can again be 
performed with a measurement in the flat region. As illustrated in Fig. 17, the 
RMS spread is still <3% in all three cases. 
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Fig. 16. For each of the shower azimuth 
angles, the peak filtered electric field 
strength per E em is constant over all en- 
ergies and particle types at a distance 
of 175 m from the shower core. 
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Fig. 17. Regardless of shower azimuth, 
the spread of the peak filtered electric 
field strength per E em in the flat region 
corresponds to only ~ 3%. 



The determination of X max values via the ratio of filtered radio amplitudes in 
the flat region and steep region also works regardless of the shower azimuth 
angle. As the calculation of this ratio removes the absolute strength of the radio 
signals, the correlation of X max with this ratio is (to first order) universal and 
independent of the shower azimuth angle. This is illustrated in Fig. 18, which 
contains all 540 showers of the three geometries in the same diagram. (The 
change of the fit parameters and RMS compared to Fig. 11, however, illustrates 
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that there are minor systematic differences between the three geometries.) 
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Fig. 18. The correiation of X max with the ratio of filtered peak amplitudes in the flat 
region and steep region is (to first order) independent of the shower azimuth angle. 
The fit parameters are a = 789.2, b = 0.3813 and c = 0.5025. The RMS deviation 
of the individual X mSLX values from the fit amounts to 25.8 gem -2 . 

4 The 45° zenith angle case 

When going from 60° zenith angle to 45° zenith angle, the qualitative be- 
haviour of the radio signal stays the same, but the lateral distance scales get 
significantly smaller. One effect which becomes, however, more important, is 
the inclusion or exclusion of energy deposited in the ground. For 45° zenith 
angle, the ground at 1400 m above sea level corresponds to a slant atmospheric 
depth of ~ 1240 gem -2 . Especially for photon-induced showers that penetrate 
very deep due to the LPM effect and have X max values of > 1200 gem -2 , a 
large fraction of the energy is then deposited in the ground. The importance 
of this aspect is illustrated by the comparison of Fig. 19 and Fig. 20. Signals 
from showers with very high values of X max are weighted down too much if E CSi \ 
is used for normalisation, i.e., the energy deposit in the ground is included. 
The normalisation with E em , on the other hand, works better, as the energy 
deposit in the ground is excluded. 

The quality of the different normalisations (with iV max , E ca \ and E em ) is further 
illustrated in Fig. 21. In contrast to the 60° zenith angle case, the normalisation 
with E ca i is now worse than that with iV max . 

If the normalisation with E em is kept, the qualitative behaviour is the same as 
in the 60° zenith angle case. This is illustrated by Fig. 22, which again shows 
the RMS spread of the peak filtered radio amplitude per E em as a function 
of lateral distance along the north (for showers coming from the south) for 
different windows of observing frequency. A clear minimum in RMS, caused 
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by the intersection of the normalised lateral radio profiles, is again present, 
but the corresponding flat region now lies at distances around 50 m, depending 
on the observing frequency. The achievable E em resolution is ~ 5%. 



When comparing the radio emission for observers at 1400 m above sea level 
(Fig. 22) and at sea level (Fig. 23), it becomes obvious that the changes in the 
45° zenith angle case are much more dramatic than for 60° zenith angle, as 
expected from the behaviour of the lateral profiles with the geometric distance 
between observer and shower maximum [7|. For observers at sea level, the 
flat region lies at ground distances between 50 and 150 m, which comprises 



a suitable range for densely equipped antenna arrays such as LOPES 11 
LOFAR H. 
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Fig. 19. For 45° zenith angle, the nor- 
malisation of filtered peak radio ampli- 
tudes with E c& \ does not work well for 
deeply penetrating gamma-ray induced 
showers caused by the LPM effect. 
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Fig. 20. Normalisation using E em works 
better than with E c& \, as the energy 
deposit in the ground is excluded. 
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Fig. 21. RMS spread of the 32-64 MHz filtered radio peak amplitudes normalised to 



E C3j i and E em , respectively, as a function of ground distance from the shower 



core in the azimuthal direction given by the shower axis for 45° zenith angle. 

Again, E em can be determined from a measurement in the flat region (Fig. 24) 
and X max can be deduced from a combined measurement in the flat region 
and a steep region (Fig. 25). The outliers at the extreme values of A max are 
again caused by the LPM effect in high-energy gamma-ray showers, being more 
prominent in the 45° zenith angle case than for 60° zenith angle showers. 
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Fig. 22. RMS spread of the filtered ra- 
dio peak amplitude func- 
tion of ground distance from the shower 
core in the azimuthal direction given 
by the shower axis for 45° zenith angle 
showers and observers at 1400 m above 
sea level. 
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Fig. 23. Same as Fig. 22 but for 
observers at sea level. 
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Fig. 24. A measurement in the flat re- 
gion at 37.5 m north directly yields the 
shower E em of 45° zenith angle show- 
ers. The fit parameter is a = 118700. 



5 Conclusions 
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Fig. 25. The ratio of peak field 
strengths in the flat region and steep re- 
gion yields the shower X max . The fit pa- 
rameters for the 45° zenith angle show- 
ers are a = 312.3, b = 0.6569 and 
c = 0.8062. The RMS spread amounts 
to ~ 19.3 gem -2 



Based on Monte Carlo simulations performed with CORSIKA and REAS2, we 
have shown that geosynchrotron radio emission from air showers of cosmic rays 
is sensitive to the primary particle energy and mass, and that this information 
content can be extracted from radio data on a shower-to-shower basis. 

For a given observing frequency and zenith angle, there is a well-defined lat- 
eral distance region in which the filtered radio peak amplitude is directly pro- 
portional to the energy deposited in the atmosphere by the electromagnetic 
cascade of an air shower, E em . Once the proportionality constant is known, 
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E em can be derived from a radio measurement in this so-called flat region. The 
RMS spread of the E em values due to shower-to-shower fluctuations amounts 
to less than 3% in case of an ideal radio signal measurement. Relating E em to 
the energy of the primary particle introduces an additional uncertainty of ~ 5 
to 7%. 

The peak field strength in the flat region can be combined with a field strength 
measured in a steep region, found at a larger lateral distance and/or higher ob- 
serving frequency. The ratio of these field strengths can then be directly related 
to the atmospheric depth of the air shower maximum, X max , on a shower-to- 
shower basis. The RMS spread of the X max values around the corresponding 
fit function due to shower-to-shower fluctuations amounts to ~ 15-20 gem -2 , 
not including any instrumental resolutions. 

Exploitation of these signal characteristics with an experimental measurement 
requires antenna spacings dense enough to allow reconstruction of the electric 
field strengths in the flat region. As the location of the flat region shifts to 
larger lateral distances for increasing zenith angles, application of this strat- 
egy is particularly interesting for inclined showers with at least 45°, better 
60°, zenith angle. (At smaller zenith angles, the flat region will shift to dis- 
tances very close to the shower core, and an extrapolation of the electric field 
strength to the core region will become necessary.) In the case of 60° zenith 
angle and for an observing frequency of 32-64 MHz, the flat region lies at a 
(ground-coordinate) distance of ~ 175 m in the azimuthal direction defined by 
the shower axis. Decreasing the observing frequency, if technically possible, in- 
creases this distance considerably. In addition, observations at lower heights 
above sea level increase the distance, especially for zenith angles smaller than 
60°. 

Assuming that the geosynchrotron model accounts for the dominant part of 
radio emission from extensive air showers, exploitation of these characteris- 
tics could thus make radio detection a powerful tool for the determination of 
the primary energy of cosmic rays and the depth of shower maximum of air 
showers, which in turn is related to the mass of the primary particle. 
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